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Abstract 

In this short note, I introduce to essential conceptual features and main building 
blocks of matrix element merging algorithms, operating on fixed order calculations 
both at leading order and next-to-leading order. The intention is purely pedagogical, 
i.e. to familiarize the reader with the essential basic concepts in a concise way, thus 
serving as an introduction to beginners and other interested readers. Unitarization is 
discussed briefly. The tutorial is highly biased towards transverse momentum ordered 
parton showers, and in particular towards merging schemes as they are implemented 
in the PYTHIA 8 general purpose Monte Carlo generator. 
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1 Hors d’oeuvre 

The present note shall provide a pedagogical and concise introduction to matrix element 
merging algorithms. It is intended to focus on the main ideas and building blocks. The 
schemes discussed are dubbed CKKW(-L), UMEPS, NL3 and UNLOPS (cited further below 
each in their respective section). With the exception of NL3, these algorithms are imple¬ 
mented in the PYTHIA 8 fflU event generator and are ready to be used. It is thus also 
natural that the discussion of merging here comes with a certain bias. For a wider overview 
over different implementations, taking into account also those available in other event gen¬ 
erators, such as e. g. MBPS |3], MENLOPS [1] and MEPS@NLO |5] in SHERPA |6], as well 
as for phenomenological implications, I strongly recommend a simultaneous study of other 
literature. 

For practical usage of the merging algorithms, e. g. for phenomenological studies, a 
consultation of the PYTHIA 8 introduction |2] and specifically the online manual [7] can be 



strongly 


The different merging algorithms are discussed on a conceptual level, the text focuses 
on what the different schemes are able to provide as well as the essential key concepts 
behind the operations. For explicit structure of the full algorithms, i. e. including the 
handling of data samples, the reader may be referred to the original publications (cited 


^Note that also the MLM merging scheme [8] is available in PYTHIA. The discussion of that scheme has 
been omitted to avoid confusion, in particular related to a different ’’order of steps” as compared to the 
algorithms discussed here. Similarly, FxFx [9] is available as a generalization to NLO. 
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further below each in their respective section). After discussing the algorithms, I close the 
note with a short discussion of the terminology merging and matching. 

The note does not contribute anything conceptually new. Its beneht lies fully in its 
pedagogical structure and clarity, together with its conciseness. I hope it will hnd its way 
to readers who will thereby become motivated to both employ as well as further improve 
QCD corrections in the held of high-energy particle physics. 

2 Motivation 

The merging of matrix elements with parton showers is the art of combining matrix element 
calculations with different jet multiplicities into a consistent combined framework and to 
consistently interface them with parton shower resummation. Parton showers provide 
a longstanding method to do all-order resummation inside of a Monte-Carlo calculation 
and can model the structure of the soft-collinear divergence of QCD. The parton shower 
generates consecutive splittings of hnal state partons by means of expressions which are 
generally called splitting functions, thus generating a cascade of an increasing number 
of hnal state partons. One of the most common choices for splitting functions are the 
DGLAP splitting kernels which are based on collinear factorization of matrix elements. 
Thus the parton shower alone can, even when being interfaced with only one single Born- 
type matrix element, give rise to higher jet multiplicities in the event data by generating 
emissions hard enough to pass the jet dehnition cut. In this way it properly accounts for 
the inclusive nature of the matrix element and resums higher-order terms in the strong 
coupling expansion of the perturbative series. Fixed-order calculations based on hard 
matrix elements on the other hand provide a more accurate description for a certain hxed 
hnal state parton mnltiplicity, in particnlar for hard emissions, i.e. in the phase space 
regions of well-separated hnal state partons. Matrix element merging methods make it 
possible to combine multiple hxed-order calculations in a consistent pictnre with the parton 
shower. The methods are designed to properly acconnt for potential donble-connting of 
hnal state conhgnrations and to provide Sndakov form factors assigning the necessary no¬ 
emission probabilities to make matrix elements exclusive in the desired way. All merging 
schemes have in common to perform a division of the emission phase space into mainly 
two distinct regions: One region where the parton shower is used and another region for 
the hard matrix elements. We refer to the cnt between these two regions as the merging 
scale. 


3 Merging scale definition 

For performing a consistent merged calcnlation, one needs to specify a certain merging scale 
dehnition, i.e. what kind of qnantity to use as a merging scale cnt. For usage with PYTHIA 8, 
is it a natural choice and advantageous to use a quantity directly related to the evolntion 
variable of the PYTHIA 8 parton shower. The PYTHIA 8 shower uses DGLAP functions, bnt 
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kinematically produces 2 —)■ 3 splittings where the third particle is involved as so-called 
recoiler (also often referred to as spectator) to ensure energy-momentum conservation. 
Since the PYTHIA 8 evolution scale is a relative pt distance of the emission, dependent on 
momenta of three particles, its concept can easily be generalized to three arbitrary partons 
of any hnal state and thus also be used on event data from a matrix element calculation. 
Already at this stage, I would like to anticipate one fundamental question that all of the 
merging algorithms will later have to address: 

’’Given a certain matrix element, how would the parton shower have generated this 
hnal state?” 

The actual importance of this question will be illustrated more thoroughly further 
below, but nota bene that this generalization of parts of the parton shower machinery 
to a matrix element calculation is apparent already here. Given an ordered set of three 
particles i,j,k out of a matrix element calculation, the corresponding PYTHIA evolution 
variable can be calculated. Minimizing this quantity over all ordered triplets of the hnal 
state will yield the merging scale dehnition pms which will be used in the remainder. Thus 
a parton level hnal state is said to be above the merging scale if 


t = min [p{i,j,k}\ where i E {hnal state partons}, 

and j G (hnal and initial state partons}, 
and k E (hnal and initial state partons} 


with the parton separation function 


Zijk{l — Zijk)Q‘lj if the radiator j is a hnal state parton, and 


Pijk ^ 


(1 ^ijk)Qij 


Qlj = {Pi+Pj)^, Zijk 

_ ^Pi{Pi+Pj+Pk) 


^i,jk 




if the radiator j is an initial state parton, and 


Qlj = -{Pi- Pjf 


. _ {Pi-Pi+Pkf 

O'fc {Pi+PkP 


( 1 ) 


( 2 ) 


is above the cut pms-, he. all jets are well separated and considered to be resolved jets, 
and will correspondingly said to be below the merging scale if f < Pms, he. at least two 
partons are in the phase space region to be populated by the parton shower. Note also 
that this merging scale dehnition can be used as a cut to regularize the hxed-order input. 


4 CKKW(-L) merging 

This section is meant to provide a schematic introduction to CKKW(-L) [TOlITT]. The 
CKKW(-L) merging schee combines several leading-order calculations and interfaces them 
with the shower. First, I will give a short overview over essential features of the CKKWL 
algorithm. Further below, some remarks will be made about the diherences between CKKW 
and Lonnblad’s implementation. 
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Since ultimatively the evaluation of observables is the main interest in any high energy 
physics Monte Carlo, we start by considering the CKKWL-based prediction for a general 
observable given by 


{O) = j #o|0(5+o,)Bo«;o + j 0{S+i,)B,wi+ j j 0{S+2j)B2W2+ (3) 


... + J ' J ^i^+Nj)^NWN 

j d(po j. . .j 0{S+nj) BnWn , 


(4) 


n=0 


where B„ denotes the differential leading order cross-section for the process of interest with 
n additional partons. These cross-sections are multiplied by weights given by 


Wr,. = 


^ofo Co ’po) ^0 /o (^0 ■Po) 

fn (^n :f^F) fn 


( 

TT ft ’Pi) _ ff (^»~ ’Pi) 

i=l 


Iln. 


X. 


1 1 Pi-1 1 Pi) I X Pni 

\i=i / Vi=i / 

fn (^n iPn) fn {^n ^Pn) 

'/n {xt.flF) X~f~ {x~,flF) 


Pms) (5) 


n 

>^n 

i=l 


OisiPi) x+_^ f+_^ {x+_^ ,Pi-l) 
«s(hi?) ^i-lfi-l(iiit-l’Pi) 


( 6 ) 


with the exception that the last Sudakov form factor Bs^^{xn, Pn, Pms) is not present for 
the highest jet multiplicity N. 

The Pi are reconstructed emission scales of successive emissions in the dehnition of the 
shower evolution variable. As has been intimated above, the merging algorithm processes 
matrix elements by trying to answer the question: “How would the parton shower have 
generated this final state?”. It will try to construct a sensible parton shower history, 
ideally starting from the 0-jet Born matrix element followed by n successive emissions with 
strictly ordered emission scales, thus being a valid event that could have equally well been 
generated by the transverse momentum ordered PYTHIA 8 shower in precisely that order. 

The CKKWL procedure then reweighs the matrix elements with parton shower weights 
that contain PDF ratios, the running of the strong coupling used in the parton shower as 
well as Sudakov form factors. The latter play a crucial role in making the events exclusive 
with respect to the merging scale. 

In an ordinary parton shower, the Sudakov form factors naturally emerge by the nature 
of the veto algorithm. This feature can also be exploited here, by running a trial shower 
on the reclustered event. Then, respective form factors can be generated by the following 
“modified” veto algorithm for each emission step: 
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Assuming an emission at a reconstructed scale pi, we need to generate a Sudakov of 
type (xj_i, pj_i, Pi). Starting from the reclustered state, perform a trial emission. 

Now distinguish between the two cases 

I. The trial emission has been generated at a higher scale than the reconstructed one. 
In that case, veto the event. 

II. The trial emission has been generated at a lower scale than the reconstructed one. 
In this case, keep the actual emission of the parton shower history and continue processing 
the event. 

Nota bene that the actual scale of the trial emission is discarded in any case. 

The “last” Sudakov factor Pn, Pwis) is generated using the same strategy by 

comparison to pus- As already mentioned, this factor is not generated in the highest jet 
multiplicity due to the beneht from keeping this type of events inclusive. Since there are 
no corresponding matrix elements for + 1 partons above the merging scale, these events 
are not double-counted and thus should not be removed in order to avoid the creation of 
a hole in the emission phase space. 

Let me mention only a few of the details and complications regarding the algorithm 
and the construction of a parton shower history and outline the solutions that were chosen. 

- Disconnected diagrams. Taking a hnal state of a matrix element and clustering its 
particles consistent with quantum numbers does not necessarily lead to a connected Feyn¬ 
man diagram, thus not yielding at all a valid emission history for the process. These 
disconnected diagrams are discarded, since connected emission histories can always be 
created. 

- Incomplete histories: The state cannot be completely clustered to the desired parton 
multiplicity and the history is incomplete. The clustering algorithm thus just leads to 
an unreasonable underlying Born. Incomplete paths are not considered if a (ordered or 
unordered) history can be generated. 

- Unordered emissions: Histories can be generated with scales that are not px ordered 
in the order of consective splittings, thus could not have been generated by the px ordered 
PYTHIA 8 parton shower in this way. Unordered histories are avoided as long as ordered 
ones can be created. Otherwise, an approximative treatment of Sudakov no-emission prob¬ 
abilities is performed. More details can be found in [TT| . 

Due to the nature of the algorithm, the total cross-section in this scheme acquires a 
dependence on the merging scale pus- This is clear, because weights are added to the 
calculation in the matrix element region being different from the weights that would have 
been obtained showering the Born matrix element onlyj§ 

The merging scheme was invented independently by CKKW and Lonnblad. A thorough 
comparison of both publications still evinces certain differences in the procedure of which 
at least the following two should be mentioned: 

—)■ CKKW apply analytic Sudakov factors to their matrix elements whereas Lonnblad 
utilizes a trial shower to generate this factor by means of the veto algorithm 

^ These effects are sometimes denoted ’’unitarity violations” in the literature since they violate the 
parton shower unitarity. 
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—)■ CKKW perform their clustering by using the longitudinally invariant kT-algorithm. 
Their procedure is thus fully deterministic whereas Lonnblad generates all possible parton 
shower histories, followed by a probabilistic selection. 


5 The concept of unitarization 


This section adresses the issue of a unitarized merging. Starting from the CKKWL algo¬ 
rithm, the necessary steps will be taken to extend the algorithm to a variant that conserves 
the parton shower unitarity. Low multiplicity states are supplemented with approximate 
higher orders from higher multiplicity states in order to retain the n-jet leading-order in¬ 
clusive cross-section. The method has been dubbed UMEPS [12] by the authors of the 
respective PYTHIA 8 implementation]^ The UMEPS event weight is given by 




fn {p^n iPn) fn (^n 

fn (^n fn i^n :t^F ) 


><n 


0‘a{Pi) TT 

Os (Mi?) 


(^i—1) Pi—1) Pi) 


(7) 


This differs from the CKKWL weight (equation Ej), as it does not contain the “last” 
Sudakov factor 115 ^^ (x„, pn, Pms)0 In UMEPS, this piece is instead generated by the higher- 
order multiplicities. The principle is motivated by the mechanism how an ordinary parton 
shower retains the inclusive cross-section. UMEPS makes this cancellation explicit: The 
emission part of the Sudakov form factor is generated mainly from the next-higher multi¬ 
plicity after the last emission has been integrated out. Since this part comes with a minus 
sign, the positive-weight contribution coming from a jet that has been emitted at a scale 
Pern is cancelled with a negative contribution where no emission between a maximal scale 
Pmax and pem has taken place. For example, events with two resolved jets can exactly 
cancel with events having one unresolved jet in the one-jet state. Note that this subtrac¬ 
tive treatment can be conveniently implemented using the parton shower history that has 
already been used for CKKWL. In the UMEPS scheme, the prediction for an observable reads 


(O) 


f d(f)o< 


1 

o 

Bi-s-o — / 62^.0 — ■ 

■ “ / Btv-s-o 



PMS 

J 

PMS 

J 

PMS 


+1 o(s+«) 




/ 


Bw-)-! 


PMS 


PMS 


+ 


unitarized merging algorithm has been suggested independently in [T3] for HERWIG++. 
"^Thus, note also wn = rejy where N denotes the highest included multiplicity. 
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0{S+N-1J 

C^{S+Nj) Bat 


BaT-I — J BAr^.Af-1 

Pms 


N ( N 

= Id^,, I..Jo(S+„j) !§„ - / b , 

n=0 ^ I 2 =n+l ^ 

V pMS 

where B^tt;^ = B„ has been used as concise notation and 


( 8 ) 



PMS 


(9) 


denotes the subtractive terms. As a subtlety, note that further integrations can be necessary 
since the merging scale measure of an event might drop below Pms due to the correct undoing 
of recoil effects, thus the additional terms B^^^ with m ^ n — 1 in equation ([8]). 

In UMEPS, the inclusive n-jet cross-section (for n < N) is not made exclusive by its own 
Sudakov resummation coming from n-jet events being processed by the shower. Instead it 
is made exclusive by the subtractive treatment of the n+ 1-jet inclusive cross-section. This 
feature can also be looked at by considering a recursive expansion of the ’’last” Sudakov 
form factor as 


^S+n i.Pm Pms) 


exp 

1 - 



) j dza^{p)Pi{p, z) 

dza^{p)Pi{p, z)As+^{pn, p) • 


( 10 ) 

( 11 ) 


The step from CKKWL to UMEPS can then also be understood as the following transition: 
Take the CKKWL-weighted expression, expand the Sudakov factor according to equation 
m The cross-section will now contain the n-jet maxtrix element times a splitting function. 
Go over to UMEPS by promoting this product to an exact (n -|- l)-jet matrix element, thus 
“undoing “ the collinear factorization. 

Loosely speaking, the philosophy of unitarization is ”we subtract what we add“: Since 
an integrated version of the terms from the higher multiplicity matrix elements is sub¬ 
tracted again, the algorithm preserves the inclusive n-jet leading-order cross-section by 
construction. 


6 NL3 

The NL3 merging scheme [TT] allows to incorporate one or several NLO calculations into 
the merging framework without inducing double-counting and still retaining resummation 
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effects. It can be considered as an NLO generalization of CKKWL. The discussion will be 
kept rather short. It is mainly shown for instructive reasons as an intermediate stage to 
the UNLOPS merging scheme which is discussed in section [7l 
The prediction for an observable in NL3 merging is given by 


{O) 


M „ „ 

E ■■■ 

f 0{S+mj) 1 +B^ 

/ L J —m,m+l 

m=0 

E Jd^J 

n=M+l 

• ■ -j 0{S+nj)BnWn 


B 


m+l^m 


( 12 ) 


where M < N NLO calculations have been added to our N tree-level calculations and 
Bm denotes an inclusive NLO calculation. The bracket notation [BmWmlmm+i Projects 
out certain powers of the as expansion, i.e. and whereas [ Cleans 

that corresponding powers of ag are explicitly removed. The aim of the algorithm is to 
provide NLO-correct results for exclusive n-jet observables while simultaneously retaining 
the parton shower resummation for all higher orders. Compared to the CKKWL prediction 
(equation Hj) , the important change is thus that M NLO calculations have been added with 
simultaneous removal of corresponding terms from the tree-level samples which would be 
double-counted otherwise. Merging schemes naturally act on exclusive cross-sections. Note 
that the inclusive NLO calculation is made exclusive by the inclusion of a phase space 
subtraction sample. 


7 UNLOPS 

UNLOPS [l5] restores the parton shower unitarity for NLO merging in a way similar to 
UMEPS in the leading-order case. It is now the NLO n-jet inclusive cross-section which 
we wish to conserve while simultaneously obtaining NLO-correct results for exclusive n-jet 
observables. The strategy is again to ’’subtract what we add“ from the lower multiplicities. 
As for NL3, we show the prediction formula for inclusive NLO input. 
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(13) 


Again, M < N NLO calculations have been added to our N tree-level calculations. 
Compared to UMEPS (equation [8]) , inclusive NLO calculations B^ have been introduced 
and, as for NL3, respective terms in the strong coupling expansion of the tree-level samples 
have been removed. Note that, if more than one NLO calculation is used, the restoration 
of unitarity requires 

I. an additional subtractive treatment with integrated emission also for the inclusive 
n-jet NLO calculations {n ^ 0) 

II. that the terms that have been removed from the n-jet tree-level calculation are 
similarly removed from their subtractive treatments in the lower-multiplicity bins (for all 
n that are supplemented with an NLO calculation) 

Higher multiplicity tree-level calculations will improve the description of exclusive n- 
jet observables while additionally the total inclusive m-jet cross-sections remain unchanged 
from the ones generated by the NLO calculations. NLO calculations with additional jets, 
in particular with multiple jets, are still not available for a lot of processes of interest. 
A common choice therefore is to supplement a single NLO 0-jet calculation with higher 
orders from tree level samples]^ 


8 A short remark on the terminology of merging and 
matching 

In particular in the older literature, a certain confusing inconsistency was present related 
to the terms merging and matching. As happens rather frequently in the etymological 
history of new words, it took time until a clear and well-dehned meaning could crystallize 
out for each of these two terms. In this section, I try to clarify them in the way they are 
used by most people these days. 

^The reader paying close attention will note that formula [T51 will in this case break down to the UMEPS 
formula [8] with the replacement Bg ^ Bg. 
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Two common standard approaches are frequently used in the implementation of NLO 
calculations: the POWHEG [16] method and MC@NLO [HI- The two approaches mainly 
differ in their subtractive treatment which is beyond the scope of this note. Formally, they 
are equivalent to the hrst order in the strong coupling expansion, differing only in the higher 
order terms. The term matching refers to the question how a (single) NLO calculation can 
be consistently interfaced with a parton shower. This is non-trivial, obviously process- 
dependent, and different for both NLO approaches. Merging on the other hand addresses 
the question how several hxed-order calculations can be combined. As is readily apparent 
from the discussion above, merging can be constructed in a generic process-independent 
way. 

Merging and matching are of course - even with these clear dehnitions - intimately 
related, and address essentially the same problems (e. g. avoiding double-counting of 
terms). Nevertheless, one should have in mind that this terminology has not been that 
clear all the time, and were thus in the past often (regarded from today’s perspective) 
confused or even used as synonyms for each other. 
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